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Abstract 

It is demonstrated that very general nonlinear dynamical systems covering 
all cases arising in practice can be brought down to rate equations of chemical 
kinetics 
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In recent years there has been great progress in the theory of classical dynamical 
systems of the form 

Fi{xi,. . . ,Xk), i = (1) 

where overdot denotes differentiation with respect to the time. Let us only recall 
the spectacular discovery of dynamical systems showing chaotic behaviour [1]. On 
the other hand, this part of the theory which deals with the universal properties 
of systems (1) can hardly be called satisfactory one. For example the Sixteenth 
Hilbert Problem (the problem of systematically counting and locating the hmit 
cycle of polynomial systems on the plane) posed in 1900 remains unsolved even for 
systems with quadratic nonlinearity [2]. 

In order to study universal properties of dynamical systems (1) one should first 
find some universal object connected with (1). An attempt in this direction was 
made by Kerner [3] who showed that very general dynamical systems (1) embracing 
everything that arises in practice can be brought down to polynomial systems 

Xi ^ Pi{xi, . . . ,xn), i^l,...,N, (2) 

where Pi arc polynomials in Xi, . . . ,Xn, by introducing suitable new variables into 
the original system (1). The aim of this letter is to demonstrate that the polynomial 
systems (2) resulting from the original nonlinear dynamical systems (1) can be fur- 
thermore reduced to systems of rate equations of chemical kinetics by introducing 
appropriate additional variables. Thus we show in the present paper that very gen- 
eral dynamical systems (1) are represented by graphs of the corresponding chemical 
reactions. 

We first recall the algorithm for construction of chemical rate equations from a 
given reaction mechanism. Let the species of the chemical reaction mechanism be 
Ai, . . . , A^. Elementary reactions are given by the stoichiometric equations 
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N N 

rririAi Yl ^riA, r = 1, . . . , M, (3) 

i=l 1=1 

where rriri, riri are the stoichiometric coefficients, rriri, G Z+ (the set of nonnega- 
tive integers) and kr is the rate constant; the formal sums from (3) are usually called 
complexes. According to the mass action law, the concentrations Xi of the species 
Ai, i = 1, . . . , N satisfy the following polynomial system 

M N 

= J2('^ri - mri)kr JJ ^7^' : i^l,...,N. (4) 

r=l j=l 

Consider now the general polynomial system (2). One finds easily that it can be 
written in the form 

UN Ji N 

n ^T' - E As n ' i^l,...,N, (5) 

where f3is G R+ (the set of positive real numbers), the vectors rij^, G Z;^ with 
coordinates riirj and rriisj, respectively are pairwise different and Ip, Jq G Z+ (it is 
understood that if Ip or = then the corresponding sum vanishes) . The criterion 
for the polynomial system (5) to have an underlying chemical reaction mechanism 
has been found by Hars and Toth [4]. Namely, they observed that whenever the 
condition 

TTiisi > if Jj > 0, s = 1,. . . ,Ji (6) 

holds for arbitrary i = 1, . . . , N, then (5) is the system of rate equations correspond- 
ing to the following sequence of chemical reactions 

N N 
^ ] f^irjAj ~^ ^ ^ fT^irjAj ~\~ A^ 

(7) 

'^isjAj > "^isjAj — Ai . 
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It can be easily checked that the general system of rate equations (4) fulfils (6). 
Therefore, (6) is the necessary condition as well. It should also be noted that there 
can exist mechanisms different from (7) having the same kinetic rate equations 
(5). Nevertheless, the advantage of (7) is that it can be constructed quickly and 
algorithmically. 

We now demonstrate that the general polynomial system (5) can be easily re- 
duced to this one satisfying the condition (6) by introducing suitable new coordi- 
nates. Consider the system (5). Suppose that the condition (6) is not valid for 
i — ii, . . . ,iq & {1, . . . , N}. On introducing the new variables of the form 

XN+p^l/xip, (8) 

we arrive at the following -|- g-dimensional system satisfying the condition (6) 

UN J, N 

Xi^Yl n ^J"' - ^Pis\{ ^7"', ^ e {1, ■ ■ ■ , \ {k, ■ ■ ■ , iq}, 

r=l j=l s=l j=l 

lip N '^ip N 

K = 1^ Q^ipfe n ^r*"'' ~ xip^N+p ^ipi n ^r^'' > (q) 
k=\ j=i 1=1 j=i 

(^ip N Jip N \ 

j2 ^ipk n ^7'"' ~ J2 (^ipi n ^7'"" ' p = 1, • • • , 
k=l 3=1 1=1 3=1 I 

We have thus shown that the general polynomial systems (5) can be brought down 
to these ones satisfying the condition (6) and the general dynamical systems (1) 
can be reduced to rate equations of chemical kinetics. An alternative method for 
the reduction of polynomial systems (2) to the chemical format has been recently 
reported by Samardzija et al [5]. Nevertheless, this method seems to be more com- 
plicated. The system (9) can be furthermore reduced to quadratic one by successive 
grouping x^Xfe, j'. A; = 1, . . . , A/" -|- 5 as a single new variable [3]. Such reduction of 
polynomial systems (5) to quadratic rate equations is the ultimate one. Indeed, if 
one persists in grouping XjX\^ as new variables then the linear infinite dimensional 
system is obtained (the nonlinear system (9) is embedded into infinite linear system 
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- such embedding is called the Carleman one [6]). The possibility of asymptotical 
reduction of general polynomial systems (2) to quadratic rate equations describing 
conservative chemical reactions without autocatalytic steps was reported by Ko- 
rzukhin [7]. Besides the perturbative nature the flaw of Korzukhin's reduction is 
its complexity. For example such simple equation as x — x is reduced (asymptoti- 
cally) within the Korzukhin approach to rate equations corresponding to chemical 
reactions of the above mentioned type by introducing five new variables. 

Example. Consider the Lorenz system 

Xi — ax2 — crxi, 

±2 = rxi — X2 — X1X2,, (10) 
is = X1X2 — bx3 . 

Owing to the term — Xix^ the system (10) is not "chemical" one. Taking into 
account (8) we introduce a new variable such that 

X4 = 1/X2. 

The system (10) implies then the following system of rate equations 

Xi — 0x2 — crxi, 

±2 = rXi - X2- XiX2X2,X4^, 

(11) 

xz = X1X2 - bxs, 

X4 — X4 + Xix^xl — rxixl . 

The sequence of chemical reactions (7) corresponding to (11) is of the form 
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Ai+A2 + A^ + A^ Ai+A^ + A^ 

Ai + A3 + 2A4 i + A3 + 3A4 
Ai + 2A4 A Ai + A4 

Here, designates the zeroth complex and the reaction A ^ (0 A) means 
that the substance A is removed from (supphed to) the chemical reactor. In order 
to reduce (11) to the system of quadratic rate equations we now introduce the 
following new variables 

X5 = xiXi, xq = (12) 

Making use of (11) one finds easily that in this variables the vector fields Xi, where 
i 7^ 4 do not depend on X4 and therefore the variable X4 can be omitted from the 
six-dimensional system obtained via (11) and (12). On renumbering — > X4 we 
finally arrive at the system of quadratic rate equations such that 

Xi — ax-i — axi, 

X2 = rxi - X2- X2X4, 

X3 = X1X2 - bxs, 

X4^ — ax3 + x1 + x1 — {a + b — l)x4 — rx^x^, 
X5 = a + x^xr, — (a — l)x5 — rx\, 

where x^ — X1X3/X2, X5 — Xi/x2- The corresponding chemical reaction mechanism 
(7) can be written as 



A1 + A2 + A3 

u 

A1 + A2 A Ai 

u u 

A2 ^ ^A3. 
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^4 + 2^5 2^5 

it 4 

^2 + ^4 A ^4 ^ >l3 ^ >l3 + >l4 

Ai + A2 + A^^Ai + A2 



2Ai A 2Ai + A4 

2^4 A 3A^ 



(13) 



In summary, wc have demonstrated in this work that very general classical dy- 
namical systems can be reduced to rate equations of chemical kinetics. The cost 
of such reduction is an increase in dimension of systems resulting from introducing 
new variables and restrictions imposed on initial data by the requirement of non- 
negativity of concentrations. In spite of these limitations, the point remains that 
general dynamical systems can be viewed as kinetic rate systems and thus they are 
represented by graphs of the corresponding chemical reactions. The problem of pos- 
sible applications of this result remains open. On the one hand, working in the field 
of nonlinear models will find that existing surprisingly strong theorems relating the 
dynamics of rate equations with the structure of the underlying graph of a chemical 
reaction mechanism such as zero defficiency theorem [8], the theorem on knots of 
trees [9] and Volpert theorems [10] (see also [11]) cannot be used in the case of 
graphs such as (13) corresponding to reaction mechanisms with autocatalytic steps. 
On the other hand, as it was mentioned by a referee, a chemist is most interested 
in inverse problem: "how to find a simple model underlying the complex reality of 
a typical chemical mechanism" (in the light of the observations presented herein it 
is really difficult job). Nevertheless, maybe at some distant stage in the future our 
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experimental ability will be such that it would be possible to "design" particular 
chemical reactor that fit given scheme and, for example, one would predict behaviour 
of a complex mechanical system from the colour of a chemical mixture. 
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